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Auto CPAP 

FIELD OF THE INVENTION 

5 The present Invention relates to method and an apparatus for the detection 
and treatment of disordered breathing during sleep, in particular to a 
method and apparatus employing an artificial neural network. 

BACKGROUND OF THE INVENTION 

10 

U.S. Patent No. 5,953,713 (Behbehani et al.) discloses a method for 
treating sleep disordered breathing comprising measuring a respiration- 
related variable at an interface placed over a patient's airway coupled to a 
pressurized gas, feeding frequency data obtained from the respiration 

15 related variable(s) into an artificial neural network trained to recognize 
patterns characterizing sleep disordered breathing; supplying pressurized 
gas to the patients airway in response to recognition of the artificial neural 
network of sleep disordered breathing. The sampling frequency of the 
pressure transducer's output disclosed in the preferred embodiment Is 512 

20 Hz. A Fourier transform is calculated every 1/16 second using a 32 sample 
values window. 

Another aspect of frequency analysis is that, on the one hand, the precision 
Is proportional to the number of input data but that, on the other hand, the 

25 response time is correspondingly increased. While high precision is welcome 
since rather small changes in breath pattern can be detected, a slower 
response increases the risk of progressive deterioration of the airway 
aperture, and thereby more severe respiratory disturbance before the 
patient is aroused. Other methods of detecting sleep disorder are based on 

30 breath-by-breath analysis 



Alternatively, if adequate treatment is not installed, the patient will be 
aroused in a more extended time perspective. 
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OBJECTS OF THE INVENTION 

It is an object of the present invention to provide an improved method for 
automatically supplying continuous positive airways pressure to a patient 

5 

It is another object of the present invention to provide an automatic 
continuous positive airways pressure apparatus (ACPAP) which lacks at 
least some of the drawbacks of state-of-the-art apparatus. 

10 Additional objects of the invention are evident from the following short 
description of the invention, the attached drawings illustrating a preferred 
embodiment, the detailed description thereof, and the appended claims. 

SUMMARY OF THE INVENTION 

15 

The present invention Is based on the insight that a direct analysis of the 
flow signal is more specific than an analysis of disordered breath, in 
particular flow limited breath, based on frequency analysis. 

20 According to the present invention is provided an automatic continuous 
positive airways pressure apparatus (ACPAP) in which the air flow from a 
CPAP or other system providing positive air pressure to a patient is 
measured for calculation of a number of parameters specific to the signal. 
The set of parameters comprises cepstrum coefficients and energy content, 

25 and is selected to indicate an apneic event of breathing during sleep, such 
as apnea, hypoapnea, and flow limitation. Data for these parameters 
collected from a large number of patients were used to train an artificial 
neural network to teach the system the variation ranges of the parameters 
for subsets of patients under a number of circumstances. The result from 

30 the artificial neural network is obtained as a low-dimensional grid of nodes 
in which each respiration type is represented by trajectory or a subsets of 
nodes. A trajectory for a normal breath looks very different from that of a 
breath during disturbed sleep. 




3 

If breathes symptomatic of a condition of disturbed sleep are detected the 
CPAP pressure is Increased. In contrast, CPAP pressure is reduced in a 
normal condition. 

5 Thus, according to the present invention is disclosed a method for the 

detection and treatment of disordered breathing during sleep employing an 
artificial neural network in which data related to breathing gas flow are 
analyzed in an artificial neural network. 

10 Specifically, the method according to the present invention comprises the 
following steps: 

• placing a mask with a tube over a patient's airway, the mask being in 
communication with a source of a pressurized breathing gas controlled 
by a CPAP, thereby establishing a respiratory circuit; 

15 • periodically sampling the gas flow in the circuit; 

• periodically calculating values for one or several parameters distinctive 
of a breathing pattern; 

• periodically feeding said parameter values to an artificial neural network 
trained to recognize breathing patterns characteristic of sleep disordered 

20 breathing; 

• analyzing said parameter values in the neural network; 

• controlling pressurized breathing gas supply in response to the output 
from said neural network. 

25 It is preferred to feed said parameter values to the network at a frequency 
of from 2 Hz to 30 Hz, preferably of about 20 Hz. 

It is preferred for said parameters to comprise cepstrum coefficients and 
energy slope. 

30 

According to a first preferred aspect of the Invention the artificial neural 
network is trained with data collected from a large number of patients. The 
data will have been collected from patients differing in many aspects: sex, 
age, body weight, breath pattern, etc. In addition, variants of sleep 
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disordered breathing such as those occurring preferentially in the back 
position, those occurring during particular stages of sleep, and those 
occurring under the influence of drugs or alcohol need to be addressed. 

5 Such data are advantageously collected in sleep laboratories in which the 
state of sleep is followed as well as the type and severity of the breathing 
disturbance is monitored by use of a polysomnography system. The 
collected data form a primary database. During the training of the artificial 
neural network the data is quantified under formation of a small secondary 
10 dedicated database which can be stored in a ACPAP. Thus, according to the 
present Invention, a dedicated secondary database obtained from a primary 
database comprising data collected from a large number of persons is 
stored in the ACPAP. 

15 According to a second aspect of the invention it is preferred to periodically 
sample the gas flow during breathing. 

The ANN comprises a number of nodes representing sets of training data. 
Each note reflects a state or an incident (feature). Neighboring nodes 

20 represent incidents of small geometric distance. In the same way as in 
training an incident vector is extracted for each flow data sample. The 
Euclidean distance from the incident vector to each node is calculated. The 
node in closest proximity to the vector is associated with it. Sequences of 
incident vectors are followed as sequences of nodes in the artificial neural 

25 network. It can be said that a sequence of nodes is the response of the 
network. Thus a trajectory in the geometric structure of the network 
(response) is followed rather than in the parameter space. The fact that the 
dimension of the network most often is smaller than the parameter space is 
of advantage since calculation thereby is simplified. The response from the 

30 network forms the basis for distinguishing between apnea, hypoapnea and 
a normal breathing state and thus, for CPAP pressure control. 

The invention thus is based on the use of an artificial neural network (ANN) 
of Kohonen-map type (associative memory; T. Kohonen, Self-Organization 
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and Associative Memory, 2 nd Ed. r Springer Vert., Berlin 1987) for detecting 
apnea or apnea-like episodes. The ANN is trained with data obtained from a 
number of patients in a sleep laboratory. The readily trained ANN forms a 
global (universal) structure of data stored in a non-volatile memory in an 
5 ACPAP. In use the breathing pattern of a patient forms trajectories (traces) 
in the ANN. A normal breathing cycle forms a dosed trajectory. A trajectory 
deviating from normal is Indicative of a breath disturbance. The ANN Is 
structured in way so as to make certain areas represent initiaf stages of 
apnea. The passage of a trajectory through such an area or touching its 
l o border indicates that the amount of air provided to the patient should be 
increased so as to re-establish normal breathing. Once breathing has been 
normalized the adduced amount of air is reduced to normal, I.e., to the pre- 
estabilshed base vatue. 

15 The artificial neural network is trained in two phases described in P. Brauer 
and P. Knagenhjelm, Infrastructure in Kohonen Maps, Proc. IEEE ICASSP, 
Glasgow 1989. 

The purpose of the analysis Is to extract, from the series of air flow rate 
20 measurements, values of the parameters chosen to classify and detect 
apneic and hypoapneic states. In each single analysis the parameters are 
made to form an incident or feature vector on which all training and 
decision-making is based. All sample values are individually analyzed in 
preparation for a quick response to changes in flow which are typical 
25 forewarnings of an apneic or hypoapneic state. 

According to a third preferred aspect of the invention linear predictive 
coding Is used to analyze the parameter values fed to the neural network. A 
linear predictive coding analysis comprising four parameters is carried out 
30 for all samples. In particular, the so-called A-parameters from the analysis 
are converted to cepstrum parameters for optimal correlation between 
parameter distance and conceptual distance, that is, so-called associativity. 
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According to a fourth preferred aspect of the invention the prediction error 
in calculating linear predictive coding is used as a basis for the parameter 
next In line. The error is filtered to counteract short-term variations and 
normalized with the total energy of the analytical window. 

5 

For calculations of energy a larger window than for the linear predictive 
coding analysis is used. The energy of the latest windows can be used to 
calculate a line the inclination which describes a trend. The difference in 
trend Is used as a further parameter. Thus, according to a fifth preferred 
10 aspect of the invention, the inclination of a trend line calculated from 
measurements and Is used as a parameter. 

According to the present invention is also disclosed an apparatus for the 
detection and treatment of disordered breathing during sleep for use with a 
15 CPAP, the apparatus including a probe for sampling breathing air flow data, 
in particular on inhalation, and an artificial neural network for analyzing, 
directly or indirectly, said data to control breathing air pressure. 

According to the present invention is also disclosed a CAPAP comprising a 
20 probe for sampling breathing air flow data, in particular on inhalation, and 
an artificial neural network for analyzing, directly or indirectly, said data to 
control breathing air pressure. 

Further variations of the present Invention are disclosed in the following 
25 detailed description of a preferred embodiment thereof illustrated in a 
drawing. 

DESCRIPTION OF THE DRAWING 

30 The Invention Is illustrated by a drawing comprising several figures, 
showing: 

Fig. 1 a block diagram regarding parameter extraction; 
Fig. 2 a feature map response in regard of example 3; 
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Fig. 3 
Fig. 4 



a feature vector; 

a feature map response In regard of example 4. 



5 DESCRIPTION OF A PREFERRED EMBODIMENT 
EXAMPLE 1. General 

PRIMARY DATA ANALYSIS. The purpose of the analysis is to extract, from 
10 the series of air flow rate measurements, values of the parameters chosen 
to classify and detect apneic and hypoapneic states. In each single analysis 
the parameters are made to form an incident or feature vector on which all 
training and decision-making is based. AH sample values are individually 
analyzed in preparation for a quick response to changes In flow which are 
15 typical forewarnings of an apneic or hypoapneic state. 

Incident vector parameters. LPC-Cepstrum. For each sample a Linear 
Predictive Coding (LPC) four-parameter analysis is carried out. The so- 
called A- parameters from the analysis are converted to cepstrum 

20 parameters for optimal correlation between parameter distance and 

conceptual distance, that is, associativity. The term cepstrum introduced by 
Bogert et al. in connection with echo time series analysis designates the 
inverse Fourier transform of the logarithm of the power spectrum of a 
signal. The transformation of a signal into its cepstrum is a homo-morphic 

25 transformation, see A.V. Oppenheim and R.W. Schafer, Discrete-Time 

Signal Processing, Prentice Hall, Englewood Cliffs, NJ, 1989. Residual. The 
error of prediction in calculating LPC is used as a basis for the following 
parameter. The error is filtered to oppose short-term variations, and is 
normalized with the total energy for the analytical window. Energy slope. 

30 For calculations of energy larger windows are used than for LPC analysis. 
The energy at the most recent windows is used to calculate a line the slope 
of which describes a trend. Difference In trend. The difference In trend Is 
used as a further parameter. 
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PARAMETERS. To detect an apneic event (i.e. central/obstructive apnea, 
hypoapnea, and flow limitations) a model is used to characterize typical 
qualities and features of the flow-signal during the event. The parameters 
of the model is chosen with the aim to be as distinct, unambiguous, and 
5 informative as possible. The set of parameters shall respond to typical 
apneic events that are readily detected by physicians. In addition to be 
sensitive to apneic events, It is Important that the parameters shall be 
insensitive to features irrelevant to the task of detecting apneic events. 

THE FEATURE VECTOR. The values of the parameters are compiled to form 
10 a vector, below named the Feature Vector. Each time the flow-signal is 
measured (sampled), the values of the Feature Vector are extracted. This 
means that if the flow signal is measured f a times per second, and N 
parameter values are needed in the model, the data rate is increased from 
/, to N f m samples per second. 

15 Prior to the extraction of parameter values, the flow signal is differentiated 
(high-pass filtered) to avoid the influence of the mean signal value. The 
mean will vary with patients and/or hardware and do not contribute in the 
classification of apneic events, and is therefore removed. Each N - 
dimensional Feature Vector can be regarded as one point in a N- 

20 dimensional signal space. 

TRAINING THE NETWORK. An Artificial Neural Network (ANN) is Iteratively 
trained to organize groups or dusters of Feature Vectors with similar 
properties. The self organizing process known as Kohonen's Self-Organizing 
Feature Map [1-2] has shown great capability of performing this task. 

25 The number of clusters is defined prior to the training and is determined by 
the required resolution of the ANN. The training is Initiated by a set of M 
clusters, randomly positioned in the N -dimensional signal space. 
The database used for training is formed by compiling the Feature Vectors 
from a large number of patients with various sleep disorders and at all 

30 stages of sleep. During the training, each Input Feature Vector is compared 
to each cluster to find the one with best resemblance to the input vector. 
This cluster is voted winner, and is adjusted towards the input vector. In 
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addition, all other clusters within a neighborhood to the winner in another 
domain, the so-called map-space are adjusted towards the input vector. 
The map-space Is usually of low dimension containing one node for each 
duster in the signal-space. The nodes are arranged in hexagonal or a 
5 square lattice, and the Euclidian distance between them defines their 
internal relation. A node's neighborhood is usually defined by a 
neighborhood function and contains the set all nodes in the beginning of the 
training whereas only a few (or none) are considered neighbors at the end. 
The further away a node is to the winner in the map-space, the less the 
10 corresponding cluster in the signal-space is adjusted towards the input 

vector. Thus all adjustments are done In the signal space, while the rules of 
adjustments are defined In the map-space. 

The training time is predetermined, and an annealing function is used to 
15 "freeze" down the system causing only small adjustments at the end of the 
training. The neighborhood function creates correlation between the signal- 
space distance and the map-space distance allowing classification to be 
performed in the (low dimensional) map-space, rather than in the more 
complicated signal-space. The method described above is known as 
20 "unsupervised learning", i.e. there is no need to use classified data in the 
training procedure described above. Classification of data into various 
apneic events is a tedious task. 

When the ANN is readily trained, the dusters will represent M features of 
25 the input flow signal including normal breathing, hypoapnea, flow- 
limitations, and apnea (provided these features are represented in the 
database used for training). The response of the ANN is proportional to the 
signal distance between the input signal and all the dusters. See figure 2. 
Often this output is of less interest in the case of classification. The output 
30 is instead used to find the node with best resemblance to a classified input, 
such as normal breathing and apneic events. This is known as the labeling 
phase in the design of the ANN. Classified Feature Vectors are presented for 
the ANN, the output is observed, and the node giving the highest output is 
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labeled with the presented feature. The actual output thereafter Is the 
label rather than the response value. 



The set of clusters are now stored In the memory of the APAP to be used In 
5 runtime mode. Patient flow-data is analyzed exactly the same way as done 
in the training phase to extract the values of the parameters used in the 
model i.e. the Feature Vector. The vector is then presented to the network 
that will produce the output label (classification) which is used by the flow- 
control logic. 

10 EXAMPLE 2. Algorithm 

DATA ACQUISITION. 

Let the flow signal be a digitized version of the analog flow-signal sampled 
at /, samples/second, giving a sequence of samples 

x„ i = 0,l,K , 

15 where x, Is short for x(i- T) I.e. the sample at time instant i-T and T-\ff s . 

PREPROCESSING. To reduce the influence of individual patient variations 
and to facilitate classification stability, the signal should pass a device to 
remove the signal mean. Any kind of steep edge high pass filter can be 
employed, thus the ideal differentiator is used for simplicity. The output 
20 from the differentiator, d, (and the input to the parameter extractor) will 
then be 

d, = x, - x,_ x where / = 0,LK , and x_, = 0 

PARAMETERS. The cepstrum coefficients have shown to well model the 
frequency content of the signal using only a few parameters (low order 
25 model). In addition, the dynamics of the cepstrum coefficients facilitate 
quantization of the parameters. Often the parameters are weighted to 
produce parameters with similar variances. The cepstrum coefficients are 
derivatives of the so called A-polynomial calculated by standard Linear 
Predictive Coding (LPC). 




As the cepstrum coefficients used do not hold information about the signal 
energy, the cepstrum will be augmented with a parameter to account for 
the long term (say 10 seconds) energy variations. This parameter shall be 
5 insensitive to the absolute level of the flow signal and only reflect the 
relative fluctuations. 

Cepstrum. T o calculate P model parameters at time k , the last Wz input 
samples are used (windowing). The sample, d k , is predicted to be 

10 • 
Thus the prediction error signal is 

M 

The task is to find the set, a , which minimizes the energy of the prediction 
error signal (i.e. finding the values of a = [a,, A ,a p ] to make J* as similar to 
15 d k as possible) over all samples within the window. The optimal solution is 
found solving the matrix equation 
Ra = -r, 

known as the Yule-Walker or normal equation, using for instance the 
Cholesky algorithm. 
20 The vector a constitutes the so called A-potynomial which is transformed 
into a set of cepstrum coefficients c using the following algorithm: 

c 0 =0 

i«t n 

The cepstrum coefficients c,,K ,c p are used as the P first coefficients in the 
Feature Vector. 

25 



• %;::; Un ri.i 
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Parameter of Energ y Trend (PET} . The energy from a number of windows 
are used to calculate a trend for the energy values. If the trend indicates 
increasing energy levels, PET Is set to zero. If the trend Is decreasing, the 
point where the trend is crossing the time-axis is calculated. A non-linear 
5 transform of this value form the PET coefficient. HBK figure 

The PET is added to the Feature Vector. The complete Feature Vector is 
thus composed according to Figure 3. 

FEATURE MAP GEOMETRY AND DEFINITIONS. Let the M map nodes be 
denoted 

10 m„ i = 0,K,M-l. 

Most often the nodes are arranged in a square (2-dimensional) grid. The 
distance between two map nodes / and J , is denoted D l f and defined as 
the squared Euclidian distance {L 2 norm) between them 
E>,j = Z. 2 (m„m / ). 

15 Let the Input Feature Vector representing sample x k be denoted y k . The 
map response in node i for feature k, , is defined as: 

5»=«p{-42XP+i), 

where the signal space distance d\ is defined as 

M 

20 and w, Is some suitable weight function. 

ANNEALING FUNCTION. The task of the annealing function is to obtain an 
equilibrium at the end of the training. The principle is that large 
adjustments are allowed in the beginning of the training whereas only small 
(or zero) adjustments are allowed at the end. How the decrease 



incorporated Is not critical. Linear, exponential, and even pulsating [4] 
decay schedules are proposed in the literature. 

INITIALIZATION. Traditionally, all data driven clustering schemes, including 
ANNS, employ random positioning of the clusters in the signal space, by 
assigning (small) random numbers to the parameters. The actual values are 
not important as long as they are not identical. The ordering of the clusters 
is also at random. 

TRAINING. The iterative algorithm adjust all clusters after each input 
Feature Vector, y k , presented. The direction of the adjustment is towards 
y k , and how much is determined partly by the annealing function, partly by 
the neighborhood function. The adjustment formulae for cluster y t at time 
instant /+1 is: 

y*('+i)=y*(0+n-('X*-y t (0), 



r*-(')=/(0**<0 

and /(/)is the annealing function and g(t)is the neighborhood function. 
Various suitable functions are discussed in [3]. 

EXAMPLE 3. Hard decision A 

Let 64 map nodes be arranged in an 8x8 square grid and numbered 0 to 63 
from the top left to the low right corner. Thus for example the map distance 
Au=l,Au=4 and D 09 =2. 

A large database is recorded containing flow-measures from several 
patients during all phases of sleep. The recordings are performed at 20 Hz 
and stored on a memory disk. The database will contain normal sleep 
breathing, flow limitations, snoring, yawning, coughing, various apneic 
events, but also mask leakage and other artifacts. 
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The database is analyzed sample for sample. The 20Hz flow-signal Is first 
passed through an ideal differentiator. A rectangular window of 180 
samples is used to form basis for extracting 4 cepstrum coefficients 
(c„K ,c 4 ) and the PET parameter. Thus the Feature Vector is a 5- 
5 dimensional vector with values extracted every 50ms. 

Samples are collected from the database in a random manner as long as 
the training proceed. The number of Iterations, T, is determined by the size 
of the database, but as a rule of thumb, 10-30 iterations per sample may 
10 be an adequate number. 

The Euclidian distance to all clusters y, i = 0,K ,M-1 are calculated and 
the duster closest to the Feature Vector is voted winner and is denoted y, . 
In calculating the distance, the following weight function is employed: 

15 w=[l 2 3 4 1] 

The neighborhood function will allow all dusters to be adjusted at all times 
(i.e. the size of neighborhood is not decreased in time), but will penalize 
dusters far away from the winner y, . 



After the training phase, the map is presented with known features such as 
normal breathing, flow limitation signals etc. For each event the response, 
S, in each node of the map is calculated (see, Fig. 2). 



20 The annealing function, /(/), follows a linear decay schedule 



/(/)=0.2-(l-*/T) f = 0,K,T 



25 



=exp(-«/*/P), 



Regions with high response for normal breathing are labeled as normal 
regions; regions reacting for flow-limitations are labeled as flow-limitation 
area and so forth. 
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It is then decided which nodes that shall represent events needed to 
respond to; i.e. an alarm signal is passed to the pressure control system. In 
Fig. 2, the map response to an input corresponding to flow limitations is 
depicted. In this case, nodes X,Y,Z will probably be labeled as a flow 
5 limitation region. 

EXAMPLE 4. Hard decision B 

This example is similar to the one described in example 1, but here a 
Hamming window of 180 samples is used instead of the rectangular. 
Furthermore the cepstrum coefficients are weighted to have approximately 
10 the same variances, whereas the PET parameter is given twice the variance 
of the cepstrum. This will give the PET parameter a little more importance 
than the rest of the parameters. The following weight function is employed: 

w=[l 1.41 1,5 4 2] 

The training is carried out exactly ass before, and the map response for a 
15 typical flow- limitation signal is shown in Fig. 4. In this case, nodes X,Y,Z will 
probably be labeled as a flow limitation region. 

EXAMPLE 5. Soft decision 

In this example, the map nodes are given a number 0, 2, 5, or 10 as labels 
to indicate the seriousness of the classification result. Thus instead of 

20 presenting an alarm signal on/off to the pressure control system, the 
number is passed on. If for Instance 0 reflects normal regions and 10 
reflects apnea, the numbers can be integrated to form an overall breathing 
status classification. If the level is very high, rapid increases in pressure is 
allowed, low levels allow for a pressure decrease, and intermediate levels 

25 result in a slow increase of the pressure. 

EXAMPLE 6. Down sampling 

If the capacity of the processor do not allow for all calculations described 
above, the LPC calculation can be decimated by a factor two or four. The 
number of samples within the analysis window must then be reduced so 




16 

that the time span of the window is still about two breathing cycles. The 
resolution of the map response will not suffer from this. 



EXAMPLE 7. Pressure regulation and runtime mode 

The pressure control system will increase the pressure one step of 0.125 
5 mbar if the ANN response is positive for 50 samples In one sequence. The 
pressure will decrease with one step of 0.125 mbar if the ANN response is 
negative for 300 samples in one sequence. The pressure will not be 
changed if the ANN response is changed during those sequences. 



In a runtime mode, a sample of the flow signal is analyzed as described 
10 above (i.e. extracting the Feature Vector), and presented to the map, now 
stored in a memory bank in the APAP-unlt. There is no need to calculate the 
exponent in the expression for map response, as the function is monotonic. 



EXAMPLE 8. Program Code 



15 // function ann_detect 

// COPYRIGHT (C) 2000 PePe Research 
// GOTHEBORGH, SWEDEN. 

#lncJude <stdlo.h> 
20 #indude "ann_detect.h" 
#lndude "apne_func.h" 

Int ann_detect( const float new_sample) { 

25 float corrf 10]; 

Int ana_dlm ■> 5; 

Int win = 180; 

float apol[5]; 

float alpha; 
30 float dvec[5]; 

int respx, respy; 

int resp_en, resp_ann; 

^ cale_corr{new_sample, win, ana_dlm+l, corr); 

/* Check levels, and slope */ 
resp_en « reg_r0(corr[0], &dvec[4]); 

lpc(corr, ana_dlm, apol, ftatpha); 

40 

a2cep(apol, ana_dlm, dvec); 
dvec[l] *= 1.41; 
dvec(2] *= 1.5; 
dvec[3] *= 4; 
45 dvec£4]*=2; 




map_resp(dvec, ana_dlm, &respx, &respy); 
resp_ann = apne_dec(respx, respy); 

if ( resp.ann == 1 || resp_en — 1 ) 
5 return 1; 

return 0; 

> 

// functions for apne detection 
10 // 

// COPYRIGHT (C) 2000 PePe Research 
// GOTHEBORGH, SWEDEN. 

# define MAP_DIM_X 8 

15 #define MAP_DIM_Y 8 

#define MAP.SIZE 64 

#define DATA_DIM 5 

#define NHEAD 256 

#deflne NMB_SAVE 100 

20 # define MaxRO 180 

#deflne MaxR0_ 2 90 

#deflne SLOPE 0.006 

#lndude <stdlo.h> 

#indude <math.h> 
25 #lnclude "map.h" 

#lndude "annpar.h" 

float distance ( const float *x, const float *y, Int d ) 

30 int I; 

float dlst, t; 

t . xt0]-y[0]; 
dlst ■ t*t; 
35 for (1=1; |<d; I++ ) < 

t = x0]-y[fj; 

dlst +o t*t; 

> 

return dlst; 

40 > 



void calc_corr( float const rec_sample, int const win, Int const ana_dim, float* corr) 

static float samp_wln[180]; 

static int first = 0; 

static float prev_sample=0; 

float dlff; 

int l J; 

/* calc diff of new and prev sample */ 
dlff = rec_sample - prev_sample; 
prev_sampte « rec_sample; 

/* Insert new sample first, move all the other */ 
for ( l=wln-l; I > 0; I- ) { 

samp_win[i] = samp_win[}-l]; 

} 

samp_wln[0] * dlff; 



for ( J - 0; J < ana_dlm; j++)< 
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corrfj] * 0; 
fbr(i-0;Kwln-J;l++)< 

corrfj] += samp_win[l]*samp_wln[J+J]; 

> 

> 

/* Add an offset to avoid near /0 operations */ 
corrfO] +=1; 



Int reg_rO(float cur_jrO) 
{ 

static float rf[MaxR0]; 
static int first = 1; 
static bit Index = 0; 
inti; 

static float prev_rt> = 0; 
float k,m, rfO, dec; 

rfO - 0.9*prev_r0 + 0.1*our_r0; 
prev_rO = rfO; 

k - 0.5*(rf0 - rftMaxRO-1]); 

m = (rfO + rltMaxRO/2] + rf[MaxR0-l]V3.; 

/♦update the saved old values */ 
for ( l=MaxR0-l; I > 0; I- )< 
rfti] - rfl>l]; 

> 

rf[0] = rfO; 

/* second rule check the regression line */ 
lf(k> -1) 

dec = 1000; 

else 

dec « -MaxR0__2*((m/k)+l); 
*out = exp(-SLOPE*dec); 

/* first rule, If energy less then 15 make alarm */ 
If(cur_r0 < 15 ){ 
return 1; 

> 

/* second rule, if regression line crosses zero line too soon make alarm * 
If C dec < 50 ){ 
return 1; 

> 

return 0; 



void 1pc(float* const acf, int const ana_dlm, float* aCoeffs, float* alpha) 

float err; 

Int I, J; 

float refls[5]; 

float sum, refl, tmp; 

err = acf[0]; 

for (i = 0; I < ana_dlm; I++) < 
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if (err <= 0) < 
reflsp] = 0; 
aGoeffs[i] - 0; 
>else{ 

sum = 0; 

for (J = 0; J < I; J++) 

sum +» aCoeffs[j] * acffl - J]; 
refl = (acffl + 1J - sum) / err; 
refls[i] = refl; 
for Q = 0; J < I / 2; J++) { 

tmp = aCoeffsrj]; 

aCbeffsrj] — refl * aCoeffc[l - J - 1]; 
aCoeffs[f - J - 1] -= refl * tmp; 

> 

if(i&i) 

aGoeffsfl / 2] -= refl * aCoeffefl / 2]; 
aCoeffs[l] a refl; 
err *■= 1 - refl * refl; 



for (I » 0; t < ana_dlm; i++) { 
aCoeffsO] = -aCoeffs[i]; 

> 

♦alpha = sqrt(err); 



void a2cep(float* const a, Int const ana_dlm, float 41 cep) 
30 < 

Int I, J, Ind; 
float sum; 

cep[0] = -l*a[0]; 
35 for (I = 2; I < ana_dlm+l; I++) { 

sum - l*a[l-l]; 
for(J-2;J<=l;j++)< 
ind = 

sum » sum + arj-2]*cep[lnd-l]*ind; 

40 } 

cep[l-l] = -sum/I; 

> 



void map_resp(float* const dvec, Int const ana_dlm, Int* respx, Int* respy) 
< 

int xm, ym, Ind, I; 
float disc; 

float min_dist = 1000000; 
for ( xm - 0; xm < MAP_DIM_X; xm ++ ) { 
for ( ym * 0; ym < MAP_DIM_Y; ym ++ > { 

Ind = (xm + (ym*MAP_DIM_X))*DATA_DIM; 
if ( (dlst « distance(dvec, &map[lnd], DATA_DIM)) < mln_dist ){ 
* respx = xm; 
♦respy = ym; 
mln.dlst m dlst; 

> 



> 
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Int apne_dec< Int const respx, bit const raspy, float* d_out) 
{ 

static float mean_x « MAP_DIM_JV2; 
5 static float mean_y o MAP_DIM_Y/2; 

static int first = 0; 
static float resp_hist[NMB_SAVE]; 
static int n_hlst_x = 0; 
static Int n_hlstL_y = 0; 
10 float racLfllt; 

float rad_lnst; 
double dx, dy, dfx, dry; 
int current_ln_area « 0; 
Int new_ln » 0, 1; 

15 

dx = (respx-area_jc_mld)*(respx-area_x_mld); 
dy = (respy-area_y_mkj)*(respy-area_y_mld); 
radjnst = sqrt(dx + dy); 

20 mean_x - lp_a*mean_x + (l-lp_a)*respx; 

mean_y = lp_a*mean_y + (l-lp_a)*respy; 

dfx « (mean_x-x_cen)*(mean_x-x_cen); 
dfy = (mean_y-y_cen)*(mean_y-y_cen); 
25 rad_fllt = sqrt(dfx +dfy); 

/♦check If the current response Is In the area */ 
If ( area_rad > radjnst ) new_In = 1; 

30 for ( I m NMB_SAVE-1; I > 0; i~ )< 

resp_hlst[l] = resp_hlst[l-l]; 

> 

resp_hlst[0] - newjn; 

35 for ( I =. 0; l< nb_resp_hist; I++ ) < 

if { 1 == resp_hist[l] ) currentjn_area++; 

> 

If ( mean__x > det_x_mln && 
mean_y > det_y_mln && 
40 mean _x < det_J?_max ft& 
mean_y < det_y_max) 

/* standard case, when detecting a simple square area of the fittred response */ 
return 0; 
else If (rad > rad_fllt) 
45 /* within drde calc, otherwise same as above */ 

return 0; 

else if (resp_ln_area < current_tn_area ) 

/'instant resp totaly over the given number */ 
return 0; 

50 else 

return 1; 

>; 

55 

EXAMPLE 9. Patient treatment 



Equipment . Breas CPAP pvlO prototype (Breas AB, MOIndal, Sweden) with 
internal flow measurement; PSG system, EMBLA™ Polysomnography 
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(Embla hf, Reykjavik, Iceland); PC with artificial neural network software 
for sleep disturbance detection. 

Patients . 6 males aged from 31 to 60 years, suffering from sleep disorders 
of various kind. 

5 Measurement. The patient was set up with all sensors from the PSG 

system: EEG, EOG, EMG, oxygen saturation, pulse, nasai flow thermistor 
measurement, body position. With this setup it is possible to determine 
sleep stages, arousals, sleep apnea, sleep hypoapnea, and other sleep 
related events. The CPAP provides air through a tube via a nasal mask to 

10 the patient. Via the CPAP communication interface the PC was connected to 
the CPAP. The PC software can read air flow values from the CPAP and set 
new pressure set points on the CPAP by the communication interface. The 
information was exchanged at a rate of about 20 Hz. The PC program feeds 
the flow values into the artificial neural network (ANN). The output from the 

15 ANN is entered into a pressure regulation algorithm (PRA). The pressure 
regulation algorithm calculates a new pressure set point and activates the 
new value in the CPAP. The output from the artificial neural network and 
the pressure setpoint Is read by the PSG system. 

Evaluation . In a normal CPAP titration (sleep disorder analysis) a CPAP is 
20 connected to a patient during sleep. The CPAP pressure is adjusted during 
the night so as to put the patient in a state with no indications of sleep 
disorders. This pressure is the one used in the CPAP treatment. The 
patient's need for various CPAP pressures can be seen with the PSG system 
breath-by- breath. The required CPAP pressure varies depending on sleep 
25 stage, body position, etc.. The data from the auto CPAP test was evaluated 
in the PSG system in the same manner as for CPAP titration by a physician 
used to evaluate patients receiving CPAP treatment. Thereby the 
correlation between the detection of sleep disorder by the ANN and the 
analysis in the PSG system could be determined. The good correlation 
30 obtained indicated that the AFN reacted correctly. 



22 



References 

[1] T. Kohonen, T. Torkkola, M. Shozakal, O. Venta "Phonetic typewriter 

for Finnish and Japanese", Proc IEEE ICASSP, New York, New York, 1988 
[2] H.C. Leung, V.W. Zue "Some Phonetic Recognition Experiments using 

Artificial Neural Nets". , Proc IEEE ICASSP, New York, New York, 1988 
[3] P.Brauer "Infrastructure In Kohonen maps", Proc. IEEE ICASSP, 

Glasgow, Scotland, 1989. 
[4] P. Knagenhjelm "A Recursive Design Method for Robust Vector 

Quantization", Proc. ICSPAT, Boston, Massachusetts, 1992. 




23 

CLAIMS 

1. A method for the detection and treatment of disordered breathing during 
sleep employing an artificial neural network (ANN) in which data related to 
5 breathing gas flow are analyzed, comprising: 

• placing a mask with a tube over a patient's airway, the mask being in 
communication with a source of a pressurized breathing gas controlled 
by a continuous positive airway pressure (CPAP) system, thereby 
establishing a respiratory circuit; 

10 • periodically sampling the gas flow In the circuit; 

• periodically calculating values for one or several parameters distinctive 
of various breathing patterns; 

• periodically feeding said parameter values to an ANN trained to 
recognize breathing patterns characteristic of sleep disordered 

15 breathing; 

• analyzing said parameter values in the neural network; 

• controlling pressurized breathing gas supply in response to the output 
from said ANN. 

20 2. The method of claim 1, where in said parameter values are fed to the 
network at a frequency of from 2 Hz to 30 Hz. 

3. The method of daim 2, wherein said parameter values are fed to the 
network at a frequency of about 20 Hz. 

25 

4. The method of claim 1, wherein said parameters comprise cepstrum 
coefficients, energy slope, difference in trend. 

5. The method of any of claims 1-5, wherein the ANN has been trained with 
30 data collected from a large number of patients. 



6. The method of claim 1, wherein the ANN has been trained with data 
collected from patients during a particular stage of sleep. 
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7. The method of daim 1, wherein the ANN has been trained with data 
collected from patients resting in a particular body position during sleep. 

8. The method of claim 1, wherein the ANN has been trained with data 
5 collected from patients being under influence of drugs including alcohol 

during sleep. 

9. The method of claim 1, wherein the ANN has been trained with data 
collected from patients by use of a polysomnography system. 

10 

10. The method of claim l, wherein the ANN comprises a number of nodes 
representing sets of training data. 

11. The method of claim 1, wherein the ANN is a Kohonen-map type ANN. 

15 

12. The method of daim 11, comprising a structure of data stored in a non- 
volatile memory. 

13. The method of daim 12, wherein the non-volatile memory is comprised 
20 by an auto CPAP. 

14. The method of daim 11, wherein the ANN comprises areas representing 
Initial stages of apnea. 

25 15. The method of claim 14, wherein a closed trajectory characteristic of 
the breathing pattern of the patient is made to pass through the ANN and is 
analyzed In regard of its geometric relationship to said areas representing 
initial stages of apnea, the passage of said trajectory through such area or 
touching such area being indicative of a disturbed breathing pattern and 

30 being used to control the ACPAP so as to increase the amount of air 
provided to the patient to restore normal breathing. 
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16. The method of daim 15, wherein the additional amount of air adduced 
to restore normal breathing is reduced upon a normal breathing pattern 
having been re-established. 

5 17. The method of dalm 11, comprising the use of 

linear predictive coding to analyze the parameter values fed to the ANN. 

18. The method of claim 17, wherein, in said linear predictive coding 
analysis, so-called A-parameters from the analysis are converted to 

10 cepstrum parameters for optimal correlation between parameter distance 
and conceptual distance. 

19. The method of daim 11, wherein the prediction error in calculating 
linear predictive coding is used as a basis for determining the parameter 

15 value next in line. 

20. The method of claim 19, wherein said prediction error is filtered to 
counteract short-term variations, and Is normalized with the total energy of 
the analytical window. 

20 

21. An apparatus for the detection and treatment of disordered breathing 
during sleep for use with a CPAP, the apparatus including a probe for 
sampling breathing air flow data, in particular on inhalation, and an artificial 
neural network (ANN) for analyzing, directly or indirectly, said data to 

25 control breathing air pressure. 

22. The apparatus of claim 21, wherein the ANN is a Kohonen map-type 
ANN. 

30 23. An automatic continuous positive airways pressure apparatus (ACPAP) 
comprising a probe for sampling breathing air flow data and an artifldal 
neural network (ANN) for analyzing, directly or Indirectly, said data to 
control breathing air pressure. 
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24. The apparatus of dalm 23, wherein the ANN Is a Kohonen map-type 
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Abstract 

A method for the detection and treatment of disordered breathing during 
sleep employs an artificial neural network (ANN) in which data related to 
5 breathing gas flow are analyzed. A respiratory circuit Is established by 
connecting the patient to a continuous positive airway pressure (CPAP) 
system with pressurized breathing gas supply, the gas flow in the circuit is 
periodically sampled, one or several parameters distinctive of various 
breathing patterns are periodically calculated; the parameter values are 
10 periodically fed to an ANN trained to recognize breathing patterns 
characteristic of sleep disordered breathing and are analyzed in the 
network, the CPAP pressurized breathing gas supply is controlled in 
response to the ANN output. Also disclosed is a corresponding apparatus. 

15 
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